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ABSTRACT 

We report results of a population synthesis model that follows the evolution of 
single and bi nary stars. I n this model we include the 2HeWD merger channel, 
suggested by lHan et al.l (|2002t ). for the formation of EHB stars. The physical 
parameters of the resulting EHB stars are derived from the BaSTI database, 
and are thus realistic and observationally supported. The predictions of this 
model are in good agreement with traditional population synthesis models, 
except when the spectrum of the stellar population is dominated by binary 
stars or their products, e.g., EHB stars in the UV of ETGs. We reproduce 
successfully the observed CMD and SED of the metal rich open cluster NGC 
6791. The stellar population in this cluster may be archetypal of the stellar 
population in ETGs that show the UVX phenomenon. Our models should be 
appropriate to study the UV upturn in ETGs. 

Key words: galaxies: elliptical and lenticular - ultraviolet: galaxies - stars : 
binaries - stars: horizontal branch 



1 INTRODUCTION 

Stellar population synthesis models are a useful tool to 
interpret the spectrum and colours of the light emit- 
ted by galaxies, as well as the properties of resolved 
stellar populations in the colour-magnitude diagram 
(CMD). Crucial ingredients for these models are com- 
plete sets of evolutionary tracks that describe in detail 
the time dependence of the physical properties of stars 
of different initial mass and chemical composition in 
the Hertzsprung- Russell diagram (HRD). Equally com- 
plete spectral libraries are needed to compute the spec- 
trophotometric properties of the stars at each position 
in the HRD. Weighting the stellar spectrum by the num- 
ber of stars at each of these positions at any given 
time, we obtain the integrated spectrum of the stellar 
population at this age. The number of stars of each 
mass born whenever there is a star formation event 
is given by the initial mass function (IMF), and the 
mass of gas transformed into stars as a function of time 
in each of these events is provided by the star forma- 
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(SFR). 



See, e.g ., Fioc fc Rocca-Volmerangj 
; Chariot! |2003l ), hereafter BC03. and 



Bruzual 
iMarastonl (J2005) for details. 

Despite considerable progress in recent years in 
both the quality and quantity of the ingredients avail- 
able to the population synthesis modelers, i.e., stellar 
tracks and spectra, these models still suffer from limi- 
tations in some specific regimes. In general, the prob- 
lems arise from our poor understanding of short-lived 
and/or not well characterized phases of stellar evolution, 
e.g., the thermally pulsing asymptotic giant branch (TP- 
AGB) and the extreme horizontal branch (EHB). These 
stars are relatively bright and contribute considerably 
to the total luminosity of a simple stellar population 
(SSP ). Thanks to the wo rk of P. Marigo and collabora- 
tors |Marigo et al.ll2008l ). we now have a better under- 
standing of TP-AGB stars of intermediate mass. These 
stars contribute at least 50 % of the near infrare d (NIR) 
light in a 1-2 Gyr old SSP (|Bruzual et al.ll2013l ). 

Progress has not reached as far in the other side 
of the spectrum. Stellar evolution models explain natu- 
rally the existence of canonical horizontal branch (HB) 
stars, but the formation mechanism and the evolution- 
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ary path of EHB stars is not well understood. EHB 
stars, being less luminous but hotter than canonical 
HB stars, are almost invisible in the optical range, but 
contribute significantly in the ultraviolet (UV), espe- 
cially in old stellar populations. Even though we do not 
understand well how EHB stars form (for many years 
these stars were hypothetical), the observational fact is 
that they have been obse rved in many open clusters 
ilKaluznv fc Udalskil QUi ), globula r clusters (ICatelanI 
l2009f). and n e arby g alaxies like M32 dBrown et alJl200Ch . 
See ICatelanI (|2009h for a review of open problems and 
the status of observations of EHB stars. 

It has long been thought that EHB stars are re- 
sponsible of the UV upturn phenomenon observed in 
the spectrum of early type galaxies (ETGs). This in- 
crement in the flux emitted shortward of 2000 A, also 
kn own as UV excess (UV X), was first detected in ETGs 
bv lCode fc Welch! \m$ U_ Since then, its o rigin has been 



a top ic of debate. | Carter et al.l J201lf). iBureau et all 
(|201lf ). and ISmith. Lucev fc Carter! (|2012T ) have shown 



beyond doubt that the UV upturn is present in passively 
evolving galaxies with no sign of recent star formation. 
Thus, the UV upturn must arise from evolved stars, and 
not from residual star formation. Were the UV excess 
produced by young stars, it should show signs of evo- 
lution, e.g. the temperature of the UV spectrum should 
evolve in time. This is not s upported by observ ations 
of ETGs (|Smith et alj |2012| ; ICarter et all l201lt ). It is 
thus fundamental to include properly the evolution of 
EHB stars in population synthesis models that describe 
ETGs. 

Several mechanisms have been proposed for the ori- 
gin of EHB stars. Single star evolution provides one such 
mechanism. When a low mass star ascends the red giant 
branch (RGB), its evolution is governed by the mass of 
the stellar envelope. If the mass loss rate is low, the star 
forms a deep and optically thick envelope, and, when the 
He flash occurs, the star reaches the HB. In high metalli- 
city and enhanced He abundance stars, mass loss is more 
effective and the star may loose all, or nearly all, of its 
envelope. The star then becomes bluer and fainter, and 
reaches the hottest region of the HB, the so called EHB. 
RGB stars of all metalliciti es with envelope mass < 0.05 
Mo evolve into the EHB |Dorman. Rood fc O'Connell 
1 19931). The EHB evolutionary phase is short-lived (~ 10 s 
yr) (|Pietrinferni et all 12004 ). However, the number of 
expected EHB stars rises with age. This suggests that if 
the UV upturn is related to the presence of EHB stars, 
the strength of the upturn shou ld increase with metalli- 
city and age jSrnith et alj|2012h . 

The evolution of bin ary stars leads naturally to the 
formation of EHB stars. Irian et al.l l|2002l . 120031 ') have 
shown that there are three channels by which interact- 
ing binaries may form EHB stars, referred to as hot 
subdwarf stars in their papers: ( a) the stable Roche 
overflow (RLOF) channel, which results in the forma- 



tion of hot subdwarf binaries with long orbital peri- 
ods, (b) the common envelope (CE) ejection channel, 
which results in the formation of hot subdwarf binaries 
with short orbital periods, and ( c) the merger of two 
He white dwarfs (2HeWD) to form a single hot sub- 
dwarf star. In the RLOF channel, the donor star fills 
its Roche lobe near the tip of its first giant branch as- 
cend, experiences stable mass transfer until its envelope 
is stripped off, resulting in a naked He core with a thin 
H envelope. In the CE ejection channel, the donor star 
also fills its Roche lobe near the tip of its first giant 
branch ascend, but dynamically unstable mass transfer 
leads to the formation of a CE. The ejection of this CE 
produces a naked He core with a thin H envelope. In 
both cases, the He core ignites to produce a hot sub- 
dwarf star. In the 2HeWD merger channel, a close He 
WD pair coalesces due to angular momentum loss via 
gravitational wave radiation. The merger product ig- 
nites He to become a hot subdwarf. This binary scenario 
successfully explains the main observational character- 
istics of field hot subdwarf stars: their distributions in 
the orbital-period vs. minimun-companion-mass and in 
the effective-temperature vs. surface-gravity diagrams; 
their distribution of orbital periods and mass function; 
their binary fraction, and the fraction of hot subdwarf 
binaries with WD companions; their birth rates; and 
their space density. In a stellar population, the forma- 
tion of EHB from binary stars dominates the far-UV 
part of the population's SE P at ages above 1 Gyr (see 
Figure 9 of lHan et af1l2007f ). The fraction of EHB stars 
formed through the 2HeWD merger chan nel become s 
larger than 50% after 10 Gyr (Figure 1 of iHanl^OOsI ). 
This channel produces more massive and more luminous 
EHB stars, which are the major sources of the FUV flux 
after 3.5 Gyr (Figure 7 of lHan et al1l2007h . 

Both scenarios for EHB star formation (single and 
binary) are plausible. However, the properties of the re- 
sulting EHB star depend on the formation channel, since 
the distributions o f the basic prope r ties of their progen - 
itors are different l|Han et alj|2007t ). ISmith et all (|2012l ) 
measured the UV upturn as a function of a few spectral 
indices and stellar velocity dispersion. They find that 
the (Fi/v — i) colo ur anti-corr e lates with age, contrary to 
the predictions of lHan et al.l |2007l ). and conclude that 
binary evolution as the source of the UV upturn is in- 
consistent with the observations. 



In this paper we revisit the inclusion of binary star 
evolution in population synthesis models, paying special 
attention to the predicted UV spectrum. First, we con- 
struct evolutionary t r acks f or single and binary stars us- 
ing the lHurlev et al.l |2002l ) public code. We add to this 
code the possibility t hat EHB stars f orm via the 2HeWD 
channel proposed bv lHan et al.l (|2002r i. The physical pa- 
rameters (T e ff, L) of the resulting EHB stars are ob- 
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tained from the BaSTI database!]] From our tracks we 
compute isochrones for stellar populations that include 
both single and binary stars, and build the correspond- 
ing spectral energy distribution (SED). We compare our 
models with the BC03 models, based on single star evo- 
lutionary tracks. 

The structure of the paper is as follows. In fJ5] we 
describe how binary star evolution is included in our 
population synthesis models. In Sj3] we compare the re- 
sults of this work with previous investigations. In [jjwe 
analyze recent obs ervations of the old metal rich open 
cluster NGC 6791 (|Buzzoni et alj|2012h . contrast them 
with our model predictions, and explore what this clus- 
ter can tell us about the UV upturn in ETGs. The con- 
clusions of our work are summarized in Sj5] 



2 A STELLAR POPULATION SYNTHESIS 
MODEL INCLUDING BINARY STARS 

As indicated in §1, close binary systems may be progen- 
itors of EHB stars through mass transfer or coalescence. 
We want to explore the effects of including binary stars, 
and hence EHB stars, in population synthesis models. 
We take a Monte Carlo approach to build from scratch 
a stellar population formed by both single and binary 
stars. To include the binary star population we must 
specify: (a) the mass distribution of the primary (the 
more massive) and secondary stars in the binary pairs, 
(b) the binary fraction as a function of the mass of the 
primary star, and (c) the distributions of the orbital pe- 
riod and orbit eccentricity of the binary pairs. For these 
distributions we follow as close as possible the obser- 
vations reported in the literature, as explained in the 
following subsections. 



Table 1. Binary fraction in function of spectral type 



Spectral Binary 
Type Fraction 



Reference 



O 0.72 Mason et al. f 19981 

O-B 0.65 Preibisch et al. (1999) 

B-A 0.62 ± 0.2 Patience et al. (20021 

G-K 0.58 ±0.1 Duquennov fc Mayor (19911 

M 0.49 ± 0.09 Fisher fc Marcv ( 19921 

Late M 0.26 ± 0.1 Barsi fc Reiners (20061 



parametrization: 



£(logm) oc 



exp 



(logm — logm c ) 2 



2<7 2 



if m < 1 M Q , 



^m~ 1J , if m > 1 M Q , 

where m c = 0.08 Mq, and a — 0.69. For the lower 
and upper mass cutoffs we adopt tul ~ 0.1 Mq, and 
my — 100 M , respectively. Each star of mass Mi will 
be the primary star in a pair according to the probability 
discussed in the next subsection and listed in Table [TJ 
Stars in our initial pool that are not selected as binaries 
remain as single stars. 

In order to obtain the mass M2 of the secondary 
star, it is necessary to e stablish a mass ratio distribu- 
tion. iMilone et ail l|2012h studied the properties of pho- 
tometric binaries in 59 Galactic globular clusters ob- 
served with the HST WFC/ACS, and concluded that 
the distribution of the mass ratio q = M2 /Mi is almost 
flat. Therefore, we assume that M2 follows from a uni- 
form distribution in q: 



/(<?) = 1 



with 



q < 1. 



(1) 



It should be remarked that our procedure does not alter 
the global IMF. Adding the binary and single star mass 
distributions we recover the original IMF. 



2.1 Statistics of binary stars 

2.1.1 The binary star population 

The first step in our approach is to generate the mass 
distribution of the primary and secondary stars in each 
binary pair. First, we build the distribution of the mass 
Mi of the primary star. The values of Mi are obtained 
by populating stochastically the IMF, which we choose , 
without loss of generality, to follow the IChabrierl (|2005l l 



1 Bag of Stellar Tracks and Isochrones (BaSTI) is a ro- 
bust and fast inte r face which u ses the evolutionary tracks of 
iPietrinferni et al. 

I <2004 I 20061 1 to obtain the stellar param- 
eters of HB stars in old stellar populations, inclu ding EHB 
stars. The tracks bv lPietrinferni et al.l l l2004l200rt) are com- 
puted with realistic stellar physics. The good agreement be- 
tween their theoretical predictions and observations is clear. 
We do not include their ct-enhanced stars in our database. 



2.1.2 The fraction of binaries 

In the last decades, surveys of binary systems have 
shown that stellar multiplicity is not the same for 
all sp ectral types. For example, iDucmennov fc Mavoj 
(Il991r ) find a multiplicity fraction of 0.58 ±0.1 for so- 
lar type stars, whereas i t is . 42 ± 0.09 for M dwars 
(|Fisher fc Marcvl Il992f l. lLadal (|2006l l reviewed recent 
data on stellar multiplicity and concluded that (a) most 
(~ 69 %) of the stars in the Galaxy are single, and (b) 
the binary fraction depends on the stellar spectral type 
of the primary star in such a way that this fraction in- 
creases with the mass of the primary star. Table [T] sum- 
marizes these results. 

2.1.3 Orbital parameters 

A binary pair is also characterized by its orbital period 
and orbit eccentricity. Surveys of binary systems have 
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shown that the distribution of the orbital period is not 
uniform. We assume that the orbital period of our bi- 
nary pairs follo ws the gaussian (in log P) dis tribution of 
period found by iDuquennov fc Mayor! |l99lf) : 



/(logP) = C7exp 



(logP-logP) 2 
2<x, 2 „„ n 



(2) 



where P is the period in days, logP = 4.4, and 
ciogP = 2.3. For the orbit eccentri city e, we assume 
that it follows a uniform distribution l|Zhang et al.ll2004l . 
l2005al ): 



with 



< e < 1. 



(3) 



2.2 Evolution of binary stars 

2.2.1 Evolutionary tracks for binary stars 

One important limitation to include binary stars in 
population synthesis models is the lack of evolutionary 
tracks computed ex professo for binary stars, which are 
needed to evolve in time the binary pairs described in 
H2.1.1I In the last decade there has been some impo rtant 
developments in this ar ea, e.g., iHurlev et al. (2002). The 
BSE code developed bv lHurlev et alj <|2002h provides the 
evolution of the stellar luminosity, effective temperature, 
mass, and evolutionary phase of both members of a bi- 
nary pair. The algorithm covers all evolutionary phases 
from the zero age main sequence, up to the remnant 
stage (black hole, neutron star, or white dwarf), and 
includes the AGB phase. However, the AGB phase is 
recognized only by a sharp increase in the star radius fol- 
lowing the core He burning phase. The thermally pulsing 
ABG stage is not modeled, but its main effects on the 
long term evolution of the star are taken into account 
i|Hurlev et al.ll2000h . 

The BSE code is designed to model complex bi- 
nary systems, in addition to model all aspects of single 
star evolution (SSE). Important aspects of binary inter- 
actions, such as mass transfer, mass accretion, angular 
momentum loss, supernova kicks, CE evolution, and the 
stable RLOF are modeled in the code. The occurrence 
of one or more of these processes depends on the param- 
eters describing the binary pair, defined in II sep- 
aration (or period), eccentricity, and mass ratio. These 
mechanisms acting in a binary pair, may cause a star 
of the pair to change its evolutionary path in the HRD 
with respect to the path that it would follow as a single 
star of its initial mass. The code can compute evolution- 
ary tracks for stars in the mass range from 0.1 to 100 
M e of metallicity Z from 0.0001 to 0.03. 

The orbital period and orbit eccentricity of the pair 
are input to the BSE code. In Figure Q] we show an 
example of how the values of these parameters affect 
evolutionary tracks. In the top left panel of Figure [T] we 
show the evolution of an Algol-like system. The binary 
interaction leads to the formation of an EHB star (log 



L ~ 1 — 1.5 and log T e a ~ 4.5). At the same time, the 
secondary star (in grey) increases its mass and becomes 
a BS, later on evolving to the EHB phase. It is impor- 
tant to note that the RGB tip is not reached by any of 
the components of the binary pair. In contrast, in the 
bottom left panel (single stars of the same mass as the 
pair above) both stars reach the RGB tip at log L ~ 4. 
The single stars do not go through the EHB or BS phase. 
The top middle and right panels show two binary sys- 
tems with different orbital period (P), eccentricity (e) 
and mass. Their single star counterparts are shown in 
the bottom middle and right panels. The luminosity at 
the tip of the RGB is lower if the stars evolve in a bi- 
nary system. The origin of these differences is the mass 
transfer between the stars in binary systems. A plus of 
the BSE code is that it can follow the entire evolution of 
even the most complex binary systems in a small amount 
of CPU time (near to one second per system). The code 
is thus ideal to compute tracks for large populations of 
binaries considering realistically their interactions. 



2.2.2 2HeWD merger 

IHurlev et al] (2002) assume in their BSE code that if 
2HeWD stars coalesce, the temperature becomes hot 
enough to start the triple-a process, and the merger 
transforms into a COWD. The process of this trans- 
formatio n, i.e., the EHB p hase, is not followed in detail. 
However, lHan et all (|2002T ) conclude that under specific 
conditions the 2HeWD star merg er can lead to t he for- 
mation of a single EHB (see also IWebbinkl ll984). This 
mechanism becomes very impor tant in populatio ns older 
than ~ 3 .5 Gyr (see Figu re 3 of lHan et al]|2010l and Fi- 
gure 7 of lHan et al] |2007m . 

In order to unders tand the formation of EHB stars 
lHan et al l (|2002l , |2003f > made a detailed investigation of 
the three main formation channels that can originate 
an EHB star. These channels are: the RLOF, the CE 
ejection, and the 2HeWD coalescence. 

RLOF occurs when a star in a binary pair fills its 
Roche lobe as a consequence of its evolution, e.g. dur- 
ing the RGB phase, or when the system looses angular 
moment, causing a reduction in the diameter of the or- 
bit. If mass transfer becomes unstable, e.g., if the donor 
star is a giant and the companion can not capture all 
the mass at the rate that it is transferred, this mass is 
accumulated and forms a CE surrounding both stars. It 
is then possible that this CE may be lost by the time 
the two cores spiral in to form a close system. The third 
channel occurs when 2HeWD stars coalesce, and the re- 
sulting object ignit es He in its core. T his last event is 
not followed in the IHurlev et al.1 (|2002h BSE code, nor 
included in their tracks. 

To include the 2HeWD star merger in our mo dels, 
we follow the recipe described by lHan et al] (|2002h . §4, 
which gives the conditions for He ignition in the 2HeWD 
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Figure 1. Evolutionary tracks computed with the BSE (binary) and SSE (single) stellar evolution code of lHurlev et al]|2002t 
The top panels show tracks for binary stars with different period (P), eccentricity (e) and masses. Note the differences in the 
tracks when the same star evolves as a single star (bottom panels). The track of the secondary star of the pair is shown in gray 
colour. The stellar masses are given in units of Mq and the orbital periods are in days. 



merger. The conditions for He to ignite depend on the 
mass of the more massive component, its accretion his- 
tory, and other parameters related to its the rmal struc- 
ture. I n order to determine these parameters, lHan et ail 
l|2002l ) performed a series of calculations that allowed 
them to follow in detail the evolution of the two WDs, 
the merger, and its products. An important result of 
their work is the mass distribution of the merger prod- 
ucts. This mass is restricted to a narrow range from ~ 
0.4 to ~ 0.65 Mq. This result is not sensitive to the ste- 
llar metallicity. If a 2HeWD star merger occurs in one of 
our BSE binary star tracks, we assume that the mass of 
the product is the sum of the mass of the two merging 
WDs at the moment of the coalescence. The resulting 
single star will ignite He if its mass is in the range indi- 
cated above. 

When He ignition happens in the merger prod- 
uct, we must assign an effective temperature and a 
luminosity to the star. This is done by interpolation 
in mass in a database constructed for this purpose 
from the output of the BaSTI web tool (ver. 5.0.1). 



The resulting effective temperatures are in the range 
20000 K < T c ft < 30000 K, whereas lHan et al l (|2002T l 
quote for these stars temperatures in the range of 20000 
K to 45000 K. This diff e rence come s from the fact that 
in the lPietrinferni et all 112004120061 ) tracks, t he range of 
He co re burning masses is narrower than the lHan et al.l 
(|2002T l distribution. We do not include the lowest and 
highest val ues of mass for 2 HeWD merger products in 
Figure 13 of lHan et af] i|2002T ). For the EHB stars formed 
through other channels, we use the temperature and lu- 
minosity provided by the BSE code. 

2.3 Isochrone synthesis 

The isochrone describing the position of single and bi- 
nary stars in the HRD at time t is built by interpo- 
lation in our set of BSE evolutionary tracks for the 
specific metallicity. Blue Stragglers (BS) are present 
in these isochrones. BS are stars which are currently 
burning H in the core, but their mass is larger than 
the turn off mass. They are located above and blue- 
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Figure 2. 5, 9 and 13 Gyr isochr ones for solar metall icity derived from the binar y and single sta r evolutionary tracks computed 
with our implementation of the iHurlev et al.l {2002) BSE code, assuming the IChab ricr 2003 IMF, and the binary fraction 
described in 32.1.11 The presence of EHB and BS stars is evident. 



ward of the turn off point in a CMD. This suggests 
that BS form through a mechanism that permits the 
star to stay in the main sequence (MS) despite its 
old age. There are two possible mechanisms to form 
BS. Both require that H is replenished in the stellar 
core by chemical mixing. One possibility is the coales- 
cence of two MS stars t o form a more massive MS star 
l|Sigurdsson et al.lll994T ). Other possible mechan ism is 
mass transfer in a binary system (McCrea 1964). Here 
we focus on the mass transfer model, included in the 
BSE tracks. Collisions are i mportant mostly in hig h 
stellar density environments (|Sigurdsson et al.l 1 19941 ). 
BS have been obs erved in all stellar systems: globular 



and open clusters dMilone et al.| 20121; Sigurdsspn et al.l 
ll994l ; ICenarro et alj|20ld: Ide Marchi et al.ll2006l ). dwarf 
spheroidal gal axies dMapelli et al l 2009). ultra faint 
dwarf galaxies IjOkamoto et al.ll2012T ). elliptical and spi- 



ral galaxies, and even in our own Milky W ay galaxy 
jMonachesi et al.ll201ll ; IClarkson et"afll201ll ). 



Another important aspect of the BSE isochrones is 
the presence of EHB stars, formed as described above. 
It should be noticed that HB stars in single star evo- 
lutionary tracks do not reach as hot temperatures as 
the EHB stars, even at the lowest metallicities. In Fi- 
gure [2] we show isochrones corresponding to the evolu- 
tionary tracks computed using the recipes described in 
W2.2\ EHB and BS stars are clearly seen. The existence 
of these stars is explained naturally by the evolution of 
interacting binary stars, which are thought to be present 
in star clusters and galaxies. The possibility that EHB 
stars form via the 2HeWD merger mechanism is still de- 
bated, and this hypothesis requires more theoretical and 
observational support. 
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Figure 3. Spectral evolution of an SSP of solar metallicity 
computed with our model. The age in Gyr is indicated next 
to each SED. 



2.4 Spectral energy distributions 

The integrated SED of the stellar population is obtained 
by adding the spectrum of the star at each position in 
the HRD along the isochrone, weighted by the number 
of stars at this position given by the IMF (see BC03 for 
details on the isochrone synthesis a pproach) . We use the 
BaSeL 3.1 stellar spectral library ((Westera et al.ll2002T ) 
to compute the spectrophotometric pr operties of all the 
models discussed in this paper. The IChabrierl (|2003h 
IMF is used throughout. 

Figure [3] shows the evolution in time of the SED of a 
solar metallicity SSP computed from our BSE tracks. At 
the youngest ages the UV light is dominated by short- 
lived massive MS stars, which rapidly evolve and leave 
the MS, causing a drop in the UV flux. At ages older 
than lGyr, the evolution of binary systems triggers the 
formation of EHB stars through the 2HeWD merger me- 
chanism. Even though these stars are also formed at 
younger ages through the CE and RLOF channels, its 
presence is not apparent because massive MS stars are 
about 100 times more luminous and dominate the UV 
spectrum. When the massive stars leave the MS, the 
EHB stars are responsible for the far ultraviolet emis- 
sion. 

The iHurlev et all (|2002h evolutionary tracks show 
the expected behavior with stellar metallicity, e.g., ste- 
llar lifetime increases with metallicity; the evolution of 
single low mass stars through the RGB phase determines 
the position of the star in the HB, at fixed initial stellar 
mass, lower metallicity stars evolve into hotter HB stars 
after the He flash. However, the formation of EHB stars 
through binary interactions depends only on the orbital 
parameters of the binary pair and not on the stellar me- 
tallicity. Thus, everything else being equal, the number 
of EHB stars present in a stellar population, should not 
de pend on metallicity. Th is has been reported previously 
bv lHan et all (|2002l . |2007m . Figure 0] shows the expected 




3.0 3.2 3.4 3.6 3.8 4.0 4.2 
Log X(A) 



Figure 4. SED of SSPs of different metallicities at 12 Gyr. 
The SEDs have been shifted in the vertical direction for clar- 
ity. The UV upturn short ward of 2000 A is present at all 
metallicities. 



similarity of the SED of SSPs of different metallicity at 
12 Gyr, the UV upturn is clearly seen at all metallicities. 

In Figure [6] we compare the evolution of the UV- 
optical colours predicted by our models for SSPs of diffe- 
rent metallicities. FUV and NUV refer to the GALEX 
filters of the same name, and i and r to the SDSS filters. 
The four colours shown in the figure become bluer as the 
population ages. The behavior of the FUV-NUV colour 
is quite remarkable. This colour is dominated at older 
ages by EHB stars, and becomes very blue for the four 
higher metallicities shown in the figure. In the two lower 
metallicity SSPs, the NUV flux is increased at older ages 
by HB stars, making FUV-NUV less blue than in their 
higher metallictiy counterparts. 



3 COMPARISON WITH PREVIOUS WORK 

Most of the popul ation synthesis models available 
in the literature dFioc fc Rocca-Volmerarigel 1 19971 . 
iLeitherer et al. I ll999l . BC03. iMarastonl 120051 ) are based 
on evolutionary tracks for single stars and do not follow 
the evolution of the binary stars known to be present in 
all stellar populations. Despite this fact, these models 
reproduce surprisingly well the spectra of most galaxies 
at all redshifts that have been sampled. This has been 
used as an argument in favour of the idea that binary 
stars are irrelevant in modeling the spectral evolution 
of galaxies. However, in the last decade there have been 
several important efforts to include binary star evolu- 
tion in population synthesis models. In this session we 
compare some results of our models with previous work. 

In Figure [S] we plot the spectrum of a 12 Gyr old 
stellar population of solar metallicity computed with 
three different codes. It is clear that below ~ 2000 A 
the present model shows an upturn in the flux because 
of the presence of the EHB stars, especially compared 
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Figure 6. UV-optical colour evolution for SSPs of different metallicity. Although the appearance of EHB stars is independent 
of metallicity, the presence of an extended HB in the lower metallicity populations increases the NUV flux and makes the 
FUV-NUV less blue than for the higher metallicity SSPs. 



with the model that uses the SSE tracks, which is com- 
pletely flat in this wavelength range. The BC03 model 
shows a modest increment in the UV flux because of 
the contribution of the central star of planetary nebulae 
(CSPNe) included in this model, and not in the other 
two. 

Figure [7] analogous to Figure [6] compares the 
evolution of the UV-optical colours of our models for 
Z = 0.008 and 0.02 with the BC03 models for the same 
values of Z. In the NUV-optical colours, the BC03 mod- 
els remain red as the population ages, whereas the mod- 
els with binaries become bluer by more than one magni- 
tude. In the colours including the FUV filter the BC03 
models also get bluer with age, due to the appearance 
of the CSPNe, but much less than our present models, 
dominated in this wavelength region by EHB stars. 

IZhang et all (|2004T ) and IZhang. Li fc Hani (|2005bl ) 



also use the iHurlev et al.l |2002) evolutionary tracks 
to follow binary star evolution in their models. Howe- 
ver, we adopt the default BSE values for the physi- 
cal parameters governing the RLOF and CE channels. 
The CE efficiency parameter (oce) is taken as 3.0, the 
Reimmer mass loss coefficient (rf) is assumed to be 0.5, 
and the tidal enha nce ment parameter is taken as 0.0. 
IZhang et all (|2004 ) and lZhang. Li fc Hani J2005bl ) use a 
different set of values. The mode ls by lLi fc Han I (|200Sf ) 
also use the IHurlev et al.l (|2002T ) tracks. They assume 
the binary fraction to be 50% for all spectral types, and 
examine mostly the optical range and NIR, and com- 
pute the Lick indices. We think that our approach that 
uses observationally supported distributions of orbital 
periods and binary fractions makes our models more re- 
alistic. 

The Binary Population and Spectral Syn- 
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Figure 7. Comparison of the UV-optical colour evolution for two of our SSPs with the BC03 models of the same metallicity. 
The CSPNe dominate the FUV flux in the BC03 models at late ages and make the model become bluer, but never as much as 
our present models, dominated by EHB stars in this wavelength range. 



thesis (BPASS) code llEldridge et alj |200S| ; 
lEldridee fc Stanwavll2009l : lEldridee et al.ll201lh follows 
the evol ution of fast rotating mas sive stars in binary 
systems |Eldridge fc Stanwavll20T2T) . according to their 
own evolutionary tracks, including the modeling of the 
interstellar gas surrounding the stars. Their main goal 
is to study the effects on stellar lifetime due to the 
mass loss induced by binary interactions. Our goal is 
different, since we aim at studying how the integrated 
properties of stellar populations are modified by the 
presence of interacting binarie s. Comparing our mo dels 
wit h either th e BPA SS, the IZhang et all |2005ri ), or 
the iLi fc Han I (|2008r i models would be an important 
exercise, since they study the same process governing 
different problems. 



4 NGC 6791: AN OPEN CLUSTER WITH 
UV UPTURN 

The open cluster NGC 6791 is an interesting stellar sys- 
tem which shows several particular features: (a) it lacks 



RGB stars |Kalirai et all 120071 ) . (6) the mass distribu- 
tion of its WD cooling s equence is bimodal l|Kalirai et al.l 
l2007l ; lBedin et al.ll200ct). and ( c) show s bimodality in the 
HB morphology ( Buzzoni et al.1120121 ') . All of these pecu- 
liarities are thought to be driv en by a mechanism of en- 
hanced mass loss. Particularly, iBedin et al.l l^OOct ) have 



shown that the observed fraction of submassive WDs in 
this cluster can be naturally accounted for if ~ 34% of 
the WDs in NGC 6791 are in binary systems. Addition- 
ally, from aste roseismology obser vations by the Kepler 
space mission, IStello et alj ()201 lh found evidence of the 
existence of unresolved binaries in this cluster. It is thus 
plausible that the origin of these peculiarities is related 



© 2013 RAS, MNRAS 000,[TJfl3l 



10 F. Hernandez- Perez and G. Bruzual 




1000 2000 3000 4000 5000 6000 



MA) 

Figure 5. SED of a solar metallicity stellar population at 
12 Gyr computed with three different synthesis codes. From 
top to bottom in the UV: (a) Present models, (b) BC03, and 
(c) SED corresponding to t he evolutionary tracks computed 
with the lHurlev et al.l l| 2000 f) SSE code, (a) includes and (b,c) 
ignore the evolution of binary stars. 



Table 2. Open cluster NGC 6791 



Observational parameter Reference 



Age ~ £ 


• 3 Gyr 


Brograad et al. (2012) 


[Fe/H] = 


= 0.4 ± 0.1 


Carraro et al. (1999, 2006) 


[a/Fe] = 


= solar 


Origlia et al. (2006) 


Y = 0.3 




Broeraad et al. (2012) 


m-M = 


13.51 


Brograad et al. (2011) 


E(B-V) 


= 0.09-0.18 


Carnev et al. ( 1999) 



to the evolution of a large number of binary star systems 
present in the cluster. 

The bimodality of the HB morphology in this clus- 
ter is relevant for our study. EHB and red clump (RC) 
stars co-exist in this cluster, suggesting that a fraction 
of its RGB stars were subject to enhanced mass loss, 
evolving into EHB stars, while another fraction evolved 
normally into RC stars. To explore the HB morphology 
of this cluster, we compute a synthetic model with the 
characteristics of this cluster, following the prescriptions 
of 321 The main observational parameters of the cluster 
are listed in Table [2] The metallicity of the cluster is 
Z=0.03. 

Figure [5] shows that there is very good agreement 
between the observed (left panel) and the synthetic 
(right panel) CMD of NGC 6791 in (B,B - V). The 
position and the extension of the region occupied by 
the BS are very well accounted for. The locations of the 
RC, the RGB and the EHB are also in good agreement 
in both diagrams. A remarkable feature is that in our 
model the bimodality in the HB is evident. It is impor- 




1000 2000 3000 4000 5000 6000 



m 

Figure 8. SED of NGC 6791 computed with our model. 



Table 3. Colours of the open star cluster NGC 6791 



Color 




This model 


Buzzoni et al. (2012) 


U- B 




0.602 


0.60 


B - V 




0.956 


0.97 


V - R 




0.625 


0.60 


V - I 




1.154 


1.18 


V - J 




2.115 


2.08 


FUV - 


V 


4.074 


5.22 


NUV - 


V 


4.066 


5.01 



tant to emphasize that the points located at the right 
of the MS in the synthetic CMD represent stars that 
are members of a binary system with a BS companion. 
These stars have lost a considerable amount of mass, be- 
coming fainter and cooler as they evolve through the sub 
giant and RGB phase corresponding to their lower mass. 
This suggests that it is possible that th e mechanism 
of enhanced mass loss mentioned above |Kalirai et al.l 
120071 ; iBedin et al.l 120081 : iBuzzoni et alll2012l ). could be 
driven by binary interactions. This will explain in a nat- 
ural manner why only a fraction of the stars passing 
through the RGB experience enhanced mass loss, lead- 
ing to the coexistence of a EHB and a RC normal for 
this metallicity and He abundance. 

In our simulated cluster, the fraction of EHB stars 
with respect t o the total of He burning stars is 20%. 
iKalirai et all l|2007tl argues that ~30% of the HB stars 
are hot, which is a extreme value for a high metallicity 
cluster. However. iDorman et alj (|l995h using population 
synthesis models found that for the highest UV upturn 
galaxies a fraction of 15% - 20% hot HB stars is enough 
to produce the observed colours. 

Figure [5] shows the SED corresponding to our syn- 
thetic model for NGC 6791, computed as described 
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B-V B-V 

Figure 9. Left panel: Observed CMD of NGC6791. The s mall gray dots are the d ata from the lStetson et all 1 1 120031) photometric 
catalog, th e open triangles are the 7 5 BS candidates from lAhumada &; Lapassetl l|2007h . and the open diamonds are the EHB 
stars from iKaluznv fc UdalskH lll992h . Right panel: Synthetic CMD. The dark gray line (in both panels) shows the isochrone 
computed for the cluster parameters listed in Table [2] 



in >i2.4l As expected, the large number of EHB stars 
present in the simulated cluster contribute to the UV 
flux, and the spectrum shows a UV upturn below 
~ 2500 A. Our spectrum is comparable with t he syn- 
thetic spectrum shown by iBuzzoni et alj l|2012T) . Their 
spectrum is built adding the contribution of the stars 
in the observed CMD of Figure [§] The slope and the 



flux level is similar in both spectra in the UV below 
2000 A. In the optical range, A > 4800 A, our spec- 
trum is somewhat fainter. This is not surprising be- 
cause in our model the binary interactions prevent some 
stars from reaching the RGB tip. From the synthetic 
SED we can compute the model cluster colours over the 
entire wavelength range. Table compares the colours 
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Table 4. Colours of NGC 6791 vs. ETGs 



Color 




synthetic 
NGC 6791 


observed 
ETGs 


FUV - 


V 


4.074 


4.95 


NUV - 


V 


4.066 


4.56 


FUV - 


NUV 


-0.008 


0.4993 


NUV - 


r 


4.459 


5.3465 



computed from our synthetic SED with those listed by 
iBuzzoni et all i|20T2T ). Both sets of colours are in good 
agreement. However, the NUV - V and FUV - V colours 
are bluer in our model. The bluer FUV - V colour may 
indicate a lack of RGB stars in our model. There are 
several possibilities to explain a low number of RGB 
stars, i) If the donor star in a binary pair is in the sub 
giant or RGB phase, it may loose enough mass to not be 
able to reach the tip of the RGB. This may be due to the 
adopted CE ejection efficiency parameter, ii) The initial 
orbital period distribution or binary fractions adopted 
in this paper may not be good enough for NGC 6791. 
In reality there may be more binary interactions pre- 
venting the primary star (star 1) of the pair to reach 
the RGB. We should keep in mind that CE evolution is 
the most importa nt but the least und erstood process in 
binary evolution fPodsiadlowski 200ll). 

The study of old metal rich open clusters may pro- 
vide an important clue to the study of the UV up- 
turn in ETGs, since in the clusters it is possible to 
resolve the hot component of the HB, which could 
act as a proxy to constrain th e source of the UV ex - 
cess in ETGs. Recent work by iRosenfield et al.l l|2012h 
shows that emission from EHB stars can be responsi- 
ble of the UV upturn in ETGs. In a subsequent paper 
l|Hernandez-Perez fe Bruzuall l2013!) we consider a sam- 
ple of (~ 3400) ETGs that have been detected in both 
the SDSS /DR8 and the G ALEX /GR6 surveys, and an- 
alyze the possible scenarios that can influence the for- 
mation of EHB stars in these galaxies, determining the 
variety observed in the UV spectrum of ETGs. 

In Table [4] we compare our synthetic cluster 
colours with the typical UV-optical colour of UV 
strong ETGs. FUV-V and NUV-V are taken from 
Table 1 of IBuzzoni et all (|2012l ). which is con- 
structed from a sample of ETGs from iBureau et al.l 
Hon]) and IBuzzoni fe Gonzalez-Lopezliral (|2008l ). 
The FUV-NUV and NU V - r c olours are from 
iHernandez-Perez fe Bruzuall (|2013h . This sample 
contains ~ 340 UV strong ETGs with no signs of recent 
star formation. The values listed in the table are the 
average colours for all of these galaxies. The synthetic 
cluster defines a blue limit to the observed averages. 



5 SUMMARY AND CONCLUSIONS 

We have built a population synthesis model that in- 
cludes, using a simple approach, the evolution of binary 
stars. In this mod el we include the 2HeWD merger chan- 
nel, suggested bv lHanet af1 f2002), for the formation of 
EHB stars . An important a s pect o f our model is the 
use of the IPietrinferni et all (|2004 120061 ) evolutionary 
tracks to estimate the stellar parameters of the EHB 
stars produced via the 2HeWD merger process. 

The predictions of our model are in good agree- 
ment with other well established single star models (e.g., 
BC03) , except in those instances where the spectrum of 
the stellar population is dominated by the binary stars 
or their products (e.g., EHB stars in the UV of ETGs). 

The CMD diagram predicted by our model for the 
metal rich open cluster NGC 6791 is in very good agree- 
ment with the observed CMD of this cluster. Our model 
reproduces well the position and relative number of BS 
and EHB stars seen in this cluster (Figure EH, as well as 
i ts exp ected SED (Figure(8]l computed bv lBuzzoni et al] 

torn . 

Despite the good match between model predictions 
and observations, the formation of EHB stars via the 
2HeWD merger process is still debated. The frequency 
of occurrence of this process in real stellar populations 
remains as an open question. 
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